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Abstract 



Radio frequency pulses are used in Nuclear Magnetic Resonance spectroscopy to produce unitary 
transfer of states. Pulse sequences that accomplish a desired transfer should be as short as possible in 
order to minimize the effects of relaxation, and to optimize the sensitivity of the experiments. Many 
coherence transfer experiments in NMR, involving network of coupled spins use temporary spin- 
decoupling to produce desired effective Hamiltonians. In this paper, we demonstrate that significant 
time can be saved in producing an effective Hamiltonian if spin-decoupling is avoided. We provide 
time optimal pulse sequences for producing an important class of effective Hamiltonians in three-spin 
networks. These effective Hamiltonians are useful for coherence transfer experiments in three-spin 
systems and implementation of indirect swap and A2(C/) gates in the context of NMR quantum 
computing. It is shown that computing these time optimal pulses can be reduced to geometric 
problems that involve computing sub-Riemannian geodesies. Using these geometric ideas, explicit 
expressions for the minimum time required for producing these effective Hamiltonians, transfer of 
coherence and implementation of indirect swap gates, in a 3-spin network are derived (Theorem 1 
and 2). It is demonstrated that geometric control techniques provide a systematic way of finding 
time optimal pulse sequences for transferring coherence and synthesizing unitary transformations in 
quantum networks, with considerable time savings (e.g. 42.3% for constructing indirect swap gates). 



1 Introduction 



The central theme of this paper is to compute the minimum time it takes to produce a unitary 
evolution in a network of coupled quantum systems, given that there are only certain specified 
ways we can effect the evolution. This is the problem of time optimal control of quantum systems 
^ This problem manifests itself in numerous contexts. Spectroscopic fields, like nuclear 
magnetic resonance (NMR), electron magnetic resonance and optical spectroscopy rely on a limited 
set of control variables in order to create desired unitary transformations |[ |^ . In NMR, unitary 
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transformations are used to manipulate an ensemble of nuclear spins, e.g. to transfer coherence 
between coupled spins in multidimensional NMR-experiments or to implement quantum-logic 
gates in NMR quantum computers The sequence of radio-frequency pulses that generate a 
desired unitary operator should be as short as possible in order to minimize the effects of relaxation 
or decoherence that are always present. In the context of quantum information processing, it is 
important to find the fastest way to implement quantum gates in a given quantum technology. 
Given a set of universal gates, what is the most efficient way of constructing a quantum circuit 
given that certain gates are more expensive in terms of time it takes to implement them. All these 
questions are also directly related to the question of determining the minimum time required to 
produce a unitary evolution in a quantum system. 

Recall the unitary state evolution of a quantum system is given by 

\m >- t/(t)iV'(o) >, 

where \'ip{t) > represents the systems state vector, at some time t. The unitary propagator U{t) 
evolves according to the Schroedinger's equation 

[> ^ -iH{t)U, (1) 

where H{t) is the Hamiltonian of the system. We can decompose the total Hamiltonian as 

m 

where Hd is the internal Hamiltonian of the system and corresponds to couplings or interactions in 
the system. Hj are the control Hamiltonians which can be externally effected The question we 
are interested in asking is, what is the minimum time it takes to drive this system |l| from U{0) = I 
to some desired C/f Pi 

In ^ , a general control theoretic framework for the study and design of time optimal pulse 
sequences in coherent spectroscopy was established. It was shown that the problems in the design 
of shortest pulse sequences can be reduced to questions in geometry, like computing shortest length 
paths on certain homogeneous spaces. In this paper, these geometric ideas are used to explicitly 
solve a class of problems involving control of three coupled spin i nuclei. In particular, the focus is 
on a network of coupled heteronuclear spins. We compute bounds on the minimum time required 
for transferring coherence in a three spin system and derive pulse sequences that accomplish this 
transfer. We also derive time optimal pulse sequences producing a class of effective Hamiltonians 
which are required for implementation of indirect swap and A2{U) gates in context of NMR quantum 
computing p3| |. 

The paper is organized as follows. In the following section we recapitulate the basics of product 
operator formalism used in NMR. The reader familiar with the product operator formalism may 
skip to the next section. Section |^ presents the main problem solved in this paper. In section 
^, we recapitulate the key geometric ideas required for producing time optimal pulse sequences. 
These ideas are developed in great detail in our work In section 5, we use these geometric ideas 
to compute the time optimal pulse sequences for producing a class of effective Hamiltonians in a 
network of linearly coupled heteronuclear spins. Finally these ideas are used to find pulse sequences 
for coherence-order selective in-phase coherence transfer in three spin system and synthesis of logic 
gates in NMR quantum computing. 
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2 Product Operator Basis and NMR Terminology 



The unitary evolution of n interacting spin i particles is described by an element of SU{2"), the 
special unitary group of dimension 2"'. The Lie algebra su(2") is a 4" — 1 dimensional space, identified 
with the space of traceless n x n skew-Hermitian matrices. The inner product between two skew- 
Hermitian matrix elements A and B is defined as < A,B >= tr{A^B). A orthogonal basis used for 
this space is expressed as tensor products of Pauli spin matrices (product operator basis). Recall 
the Pauli spin matrices Ix, ly, Iz defined by 




are the generators of the rotation in the two dimensional Hilbert space and basis for the Lie algebra 
of traceless skew-Hermitian matrices 5u(2). They obey the well known relations 

ly] = ih ; [ly, Iz] = ilx ; Ix] = ily (2) 



l!=l'y=l'z=\l (3) 

where 

■ 1 
1 



Notation 1 We choose an orthogonal basis {iBs} (product operator basis), for su(2") taking the 
form 

n 

=2«-i[](/to)"-, (4) 

fc=i 

a — x,y,or z and 

/fca = 1 (g) • • • (8) /a (g) 1 , (5) 

where q is an integer taking values between 1 and n, the Pauli matrix appears in the above 
equation |^ only at the k*^^ position, and 1 the two dimensional identity matrix, appears everywhere 
except at the k*^ position, a^s is 1 in g of the indices and in the remaining. Note that we must 
have q > 1 as q = corresponds to the identity matrix and is not a part of the algebra. 



Example 1 As an example for n = 2 

q = 2 



'. product basis for su(4) takes the form 

iillx, hyi Ilz, hx, hy, l2z} 
i{'2Ilxl2x, 2Iixl2y, 2Ilxl2z 
2Ilyl2xi 2/lj,/2j,, 2Ilyl2z 
^hzhx: 2Ilzl2y, ^Iizhz-} 
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Remark 1 It is very important to note that the expression I^a depends on the dimension n. For 
example, the expression for l2z for n = 2 and n = 3 is 1 , and 1 (g) (g) 1 respectively. Also 
observe that these operators are only normalized for n = 2 as 



To fix ideas, we compute one of these operators explicitly for n — 2 



1 
1 

which takes the form 

1 " 
10 
0-10 
-1 

In this paper we want to control a network of coupled heteronuclear spins. The internal Hamil- 
tonian for a network of weakly coupled spins takes the form 

Hd = 2TT^ ViUz + 27r ^ Jijhzljz- 

i ij 

Where represents Larmor frequencies for individual spins and Jij represents couplings between the 
spins. The values of the frequencies Vi and Jij depend on the particular spins being used; typically, 
— 10* — 10^ Hz while for neighboring spins Jij — 10 — 10^ Hz. Throughout this paper, we will 
assume that the Larmor frequencies of spins are well separated {\vi — Vj \ ^ In a frame rotating 

about the z axis with the spins at respective frequencies i^i, the Hamiltonian of the system takes the 
form 

Hd = 271'^ Jijiizljz- 
ij 

We can also apply external radio frequency (rf) pulses on resonance to each spin. Under the 
assumption of wide separation of larmor frequencies, the total Hamiltonian in the rotating frame 
can be approximated by 

H = 27r^^JijIizIjz + 2'lT^^{ViiIix + Vi2liy), 
ij i 

where lix and liy represent Hamiltonians that generate x and y rotations on the i^^ spin. By applica- 
tion of a resonant rf field, also called a selective pulse, we can vary vn and Vi2 and thereby perform 
selective rotations on individual spins. In this context, we use the term hard pulse if the radio- 
frequency (rf) amplitude is much larger than characteristic spin-spin couplings. Such hard pulses 
can still be spin-selective if the frequency difference between spins is larger than the rf amplitude 
(measured in frequency units) In particular, this is always the case for the heteronuclear spins 
under consideration. In many situations, it is possible to "turn off" one or more of these couplings 
Jij. This is done through standard spin decoupling techniques, for details see ^ and appendix A. 

We now present the main problem addressed in this paper. 



hz 



1 

-1 
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3 Optimal Control in Three Spin System 



Problem 1 Consider a chain of three heteronuclear spins coupled by scalar couplings ( J13 = 0). 
Furthermore assume that it is possible to selectively excite each spin (perform one qubit operations in 
context of quantum computing). The goal is to produce a desired unitary transformation U G SU (8), 
from the specified couplings and single spin operations in shortest possible time. This structure 
appears often in the NMR situation. The unitary propagator U, describing the evolution of the 
system in a suitable rotating frame is well approximated by 

6 

U = -t{Hd + Y^ UjH, )U, U(0) = I (7) 
i=i 

where 





= 27rJi2/] 


Hi 


= 27r/ia;, 


H2 


= 27r/iy, 


H3 


= 27r/2a , 


Hi 


= 27r/2y, 


H5 


= "^ttIsx , 


He 


= ■ 



The symbol J12 and J23 represents the strength of scalar couplings between spins (1,2) and (2,3) 
respectively. We will be most interested in a unitary propagator of the form 

Where the index a, /?, 7 G {x, y, z}. These propagators are hard to produce as they involve trilinear 
terms in the effective Hamiltonian. We will refer to such propagators as trilinear propagators. To 
highlight geometric ideas, here we will treat the important case of this problem when the couplings 
are both equal (J12 — J23 — J)- Without loss of any generality we assume J > 0. 



Remark 2 Please note that it suffices to compute the minimum time required to produce the 
propagators belonging to the one parameter family 

Uf = cxpi-i9 hzhzhz), e [0,4^], 

because all other propagators belonging to the set {exp(— 16* /iQ/2/3^37)|a, /9, 7 S {x, y, 2}} of trilinear 
propogators can be produced from Up in arbitrarily small time by selective hard pulses. As an 
example 

exp(-z6' hxhzhz) = exp(-i- /iy)exp(-i6' hzhzhz) exp{i-^ hy)- 

It will be shown that finding shortest pulse sequences for these propogators, constitute an essential 
step in optimal implementations of logic gates in the context of NMR quantum computing. 



Remark 3 We first compute the minimum time it takes to produce the propagator of the above type 
using spin-decoupling. The main computational tool used for this purpose is the Baker Campbell 
Hausdorff formula [BCH] g. RecaU given the generators A, B, C satisfying 

[A, S] = C , [B, C] = A, [C, A] = B. 



5 



The BCH implies 

exp(At)Sexp(-At) = B cost + C sint, 

and therefore 

exp{At) exp(_B) exp{—At) — exp{B cost + C sint). 
This can be then used in problem |^ to produce a propagator of the form exp(— i0 Iizl2zlzz)- 

The standard procedure uses decoupling and operates by first decoupling spin 3 from the network 
(this can be achieved by standard refocusing techniques [Q, see Fig. |^(A). A brief review of the 
basic ideas involved in spin-decoupling is presented from a control viewpoint in appendix A). The 
effective Hamiltonian then takes the form 

Hlff = 2'kJIizI2z- 

Now by use of external rf pulses and the Hamiltonian H^jj: , we can generate the unitary propagator 
exp{—iTrIizl2x) as follows. 



exp(-i-/2y)exp(-z-^-^)exp(i-/2y) = cxpi-inlizhx)- 



The creation of this propogator takes ^ units of time 



Similarly by decoupling spin 1 from the network, we are left with an effective Hamiltonian H^ff 
2ttJI2zI3z: which can be used along with external rf pulses to produce a propagator exp(— 
which takes another units of time. Now using the commutation relations 

[2/12/2^, 2/2y/3z] = i^Ilzhzhz, 
Wlzhzhz,2hzhx\ = l^hyhz. 

Wizhzhz^^hzhx] = i^hyhz- 



We obtain that 



exp(-i7r/iz/2j:)exp(-i^^^^^)exp(i7r/i2/2a;) = eyip{~i6Iizhzhz) 



2 

Therefore the total time required to produce the unitary propagator is 

1 e 1 A-K + e 2 + K 



2J AttJ 2J 4ttJ 2J ' 

where n = ^ (see Fig ^ . 

We will show that this propagator can be produced in a significantly shorter time using pulse se- 
quences derived using ideas from results in geometrical control theory. Before we turn to time 
optimal pulse sequences, we give new implementations of the trilinear propagators that are con- 
siderably shorter than the ones given in remark |^, even though they are not time optimal. These 
sequences do not involve decoupling. We present one such sequence here, for comparison with the 
time optimal pulse sequences in theorem |^(see Fig. ||(B)). 

Notation 2 Let A = -i{Iizhx + hxhz), B = -i{Iizhy + hyhz), C = -i{2Iizl2zhz + ^) and 
D = —i{4:Iizl2zl3z)- Then observe the following commutation relations hold 

[A,B]^C; [B,C]=A- [C,A]=B. (8) 

[A,D]^~B; [B,D]^A. 
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Figure 1: The graph shows the comparison of time required by pulse sequences for creating trihnear 
propogators as a function of k = 0/27r. (a) Pulse sequence using spin-decoupling, (b) improved 
sequence without decoupling (see remark 4), (c) time optimal pulse sequence (see theorem]^). 



Definition 1 Any set of three generators A,B,C satisfying the equation (g) will be referred to as 
the so (3) Lie algebra. 



Remark 4 Using the commutation relations stated above, it follows from BCH that 

P = exp(|A) exp(^B) exp(^|A) = exp(-z0(/i,/2./3. + ^)). 

It takes arbitrarily small time to generate the propagator Q = exp(z6'-^)), using selective hard 
pulses. Thus the time required to generate the desired propagator PQ — exjp{—i{9Iizl2zl3z)) is just 
the time needed to produce P, which can be computed explicitly. The propagator exp(^A) requires 
jj units of time, and the propagator exp(|i3) requires -^-j units of time. Hence the total time is 

1 9 1 _ 1 + K 
47 ^ M ^ 47 ~ 2J ' 

Thus we see that it is possible to reduce the time of pulse sequences for implementing desired effective 
Hamiltonians, by not decoupling spins in the network. The savings are as much as 50% for small k 
(see figure 0) 



We now state results on time optimal pulse sequences for coherence transfer and synthesis of logic 
gates in 3 spin systems. The main theorems of this paper are stated as follows. 



Theorem 1 Given the spin system in (^, with J12 = J23 = J and J13 = 0, the minimmn time 
t*{UF) required to produce a propagator of the form Up ~ exp{—i9Iizl2zl3z), & G [0,47r] is given 

by 

^ v/27rg-(g/2)2 ^ ^k{A-^) 

^ ^' 2-nJ 2J 

where k = 
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This theorem can be used to compute the minimum time and the shortest pulse sequence required 
for in-phase coherence transfer in the three spin network given by equation (|^) and construction of 
swap gates between spin 1 and 3. This is stated in the following theorem. 

Theorem 2 (Indirect Swap Gates and Coherence Transfer:) Given the spin system in (|^), 
with Ji2 = J23 = J and J13 — 0, the minimum time required for producing a swap gate between spin 
1 and 3 is The minimum time required for the complete in-phase transfer — [Ii^ — iliy) to 



Remark 5 The conventional approach for the above indirect swap gate involves three direct swap 
operations. The first operation swaps spin 1 and 2, followed by a swap 2 and 3 and finally a swap 
between 1 and 2 again. Each operation takes ^ units of time. The total time for this pulse sequence 
is jj. Compared to this the time optimal sequence only takes = 57.7% of the total time. It is 

possible to transfer /f completely using two sequential selective isotropic steps that involves 

decoupling, each of which takes units of time This takes in total j units of time. The 

improved pulse sequence takes at most = 86.6% of this time. 

We now derive the time optimal pulse sequences that give the shortest times described in above 
theorems. We begin by recapitulating the main geometric ideas developed in for finding these 
time optimal pulse sequences. 



4 Main Ideas 



Let G denote the unitary group under consideration. In the equation 

m 

i=i 

the set of all U' E G that can be reached from Identity / within time t will be denoted by R(/, t). 
We define 



t*{UF) = inf {< > 0| C/f e R(/,t)} 



where R(/, i) is the closure of the set R(/,i), and / is the identity element. <*(C/f) is called the 
infimizing time for producing the propagator Up- Observe that the control Hamiltonians {-ffj}, 
generate a subgroup K, given by 

K^eM{H,}LA), 

where {Hj^LA is the Lie algebra generated by {—iHi, — ii?2, ■ • ■ , —iHm}- It is assumed that the 
strength of the control Hamiltonians can be made arbitrary large. This is a good approximation to 
the case when the strength of external Hamiltonians can be made large compared to the internal 
couplings represented by Hd- Under these assumptions the search for time optimal control laws can 
be reduced to finding constrained shortest length paths in the space G/K. It can be shown that 
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Theorem 3 (Equivalence theorem): The infimizing time t*{UF) for steering the system 

m 

from C/(0) = / to Up is the same as the minimum time required for steering the adjoint system 

P = HP, He AdKi-iHd), PeG (9) 
from P(0) = / to KUf, where Adxi-iHa) = {kl{-iHd)ki\ki e K}. 

We will use this result to find time optimal pulse sequences for 3-spin system. The key observation 
leading to the Equivalence theorem is summarized as follows. 

(Minimum time to go between cosets:) If the strength of the control Hamiltonians can be 
made very large, then starting from identity propagator, any unitary propagator belonging to K can 
be produced in arbitrarily small time. This notion of arbitrarily small time is made rigorous using 
the concept of infimizing time as defined earlier. Therefore if Up G K then t*(UF) = 0. Similarly, 
starting from Ui, any kUi, k € K can be reached in arbitrarily small time. This strongly suggests 
that to find the time optimal controls Vi which drive the evolution (|l|) from Ui to U2 in minimum 
possible time, we should look for the fastest way to get from the coset KUi to KU2 (the coset KUi 
denotes the set {kUi\k e K}). 

(ControlUng the direction of flow in G/K space:) The problem of finding the fastest way 
to get between points in G reduces to finding the fastest way to get between corresponding points 
(cosets) in G/K space. Let g represent the Lie algebra of the generators of G and i — {HjJla 
represent the Lie algebra of the generators of the subgroup K. Consider the decomposition g = p©6 
such that p is orthogonal to 6 and represents all possible directions in the G/K space. The flow 
in the group G, is governed by the evolution equation (^ and therefore constraints the accessible 
directions in the G/K space. The directly accessible directions in G/K, are represented by the 
set Adxi—iHd). To see this, observe that the control Hamiltonians do not generate any motion in 
G/K space as they only produce motion inside a coset. Therefore all the motion in G/K space is 
generated by the drift Hamiltonian Hd- Let ki and ^2 belong to K, the coset containing identity. 
Under the drift Hamiltonian Hd-, these propagators after time St, will evolve to exp{—iHd St)ki and 
exp{—iHd 6t)k2, respectively. Note 

exp{—iHd 5t)ki ~ ki{k\ exp(— ii/^ 5t)ki) 

and thus is an element of the coset represented by 

fc| exp(— ii/d 5t)ki = exp{~ik\Hdki St). 

Similarly exp{'-iHdSt)k2 belongs to the coset represented by element exp{—ik2Hdk2 6t) . Thus in 
G/K, we can choose to move in directions given by k\{~iHd)ki or kl(~iHd)k2, depending on the 
initial point ki or fc2. Therefore all directions AdK{—iHd) in G/K can be generated by the choice 
of the initial k e K , hy use of control Hamiltonians {Hj} (We can move in K so fast that the 
system hardly evolves under Hd in that time). The set AdK{—iHd) is called the adjoint orbit of 
—iHd under the action of the subgroup K. This form of direction control has been defined as an 
adjoint control system Observe that the rate of movement in the G/K space is always constant 
because all elements of AdK{iHd) have the same norm, \\Hd\\ = \\k^Hdk\\ {k is unitary so kk^ is 
identity). Therefore the problem of finding the fastest way to get between two points in the space 
G/K reduces to finding the shortest path between those two points under the constraint that the 
tangent direction of the path must always belong to the set Adfc{~iHd)- This is the content of 
equivalence theorem. 
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(Finding Sub-Riemannian Geodesies in Homogeneous spaces:) The set of accessible direc- 
tions AdK{~iHd), in general case is not the whole of p, the set of all possible directions in G/K. 
Therefore all the directions in G/K space are not directly accessible. However, motion in all di- 
rections in G/K space may be achieved by a back and forth motion in directions we can directly 
access. This is the usual idea of generating new directions of motion by using non-commuting 
generators ( exp(e^) exp(ei3) exp(— ev4) exp(— eS) ~ exp(— e^[A, i?]) ). The problems of this na- 
ture, where one is required to compute the shortest paths between points on a manifold subject to 
the constraint that the tangent to the path always belong to a subset of all permissible directions 
have been well studied under sub-Riemannian geometry. These contrained geodesies are called the 
sub-Riemannian geodesies The problem of finding time optimal control laws, then reduces to 
finding sub-Riemannian geodesies in the space G/K, where the set of accessible directions is the set 
AdKi-iHd). 

In H], these sub-Riemannian geodesies were computed for the space su{2)^^'su{2) ' '^^ ^^"^ context of 
optimal control of coupled 2-spin systems. It was shown that the space su{2y»s'ui2) struc- 
ture of a Riemannian symmetric space which facilitates explicit computation of these constrained 
geodesies. In the following sections we will study these sub-Riemannian geodesies to compute the 
time optimal control for three spin systems. 



5 Time Optimal Pulse Sequences 

In the following lemma, we describe the infimizing time for the heteronuclear three spin system, 
described by the equation with J12 = J23 = J and J13 — 0, in terms of its associated adjoint 
control system 

P^HP, He AdK{-i2TTJ{hJ2z + hzhz)), 
where K denotes the subgroup generated by control Hamiltonians {-ffjjj^i. 

Lemma 1 In equation let K denote the subgroup generated by control Hamiltonians {Hj}^^^. 
The infimizing time t* {Up), required to produce a unitary propagator Up is the same as the minimum 
time T, required to steer the adjoint control system 

P = HP, HeAdK{-i2^J{hzhz + hzhz)). (10) 
from P(0) = / to P{T) e KUp- 

Proof: The lemma follows directly from the equivalence theorem ^ Q.E.D. 

In the following theorem, we develop a characterization of time optimal control laws for the 
adjoint control system @. This characterization is obtained using the maximum principle of Pon- 
tryagin. We briefly review the maximum principle here. The reader is advised to look at the reference 
[pl for more details. 



Remark 6 Pontryagin MsLximum Principle: Consider the control problem of minimizing the 
time required to steer the control system 

X = J{x, u), xeM", uenc R'', 

from some initial state x(0) = Xq to some final state xi. The Pontryagin maximum principle states 
that if the control u{t) and the corresponding trajectory x{t) are time optimal then there exists an 
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absolutely continuous vector A(<) E M", such that the Hamiltonian function Sj{x{t), X{t),u{t)) — 
{t)f{x{t),u{t)), satisfies 

and 

A,W = -^, jel...n. 

The vector A(^) is called the adjoint vector and any triple (x, A, u) that satisfies the above conditions 
is called an extremal pair. The basic ideas of this theorem can be then generalized to control 
problems defined on Lie Groups |l2j. We use these ideas to give the necessary conditions for the 
time optimal control laws for the adjoint control system (^). 



Theorem 4 For the adjoint control system (^), if 7Y(t) is the time-optimal control law, and P{t) 
is the corresponding optimal trajectory, such that P{0) — I and P{T) £ KUf, then for t E [0,T], 
there exists M{t) e p, (directions in G/K space) such that 

Hit) = argmax ir(7i:M(t)), H E AdK{-iHa), (11) 
H 

'^^^^^ = n{t)p{t), (12) 



dt 
dM{t) 
dt 



m).M{t)] (13) 



Proof: First note TO = —Ti, as TL is skew-Hermitian. We represent the linear functional on P as 
4>\{P) = tr{X^T-CP) with PX^ E p (the directions corresponding to G/K space). The Hamiltonian 
function is then 

^{Pit), x{t), n{t)) = tr(At {t)n{t)P(t)). 

Then the maximum principle gives 

n{t) = argmax triHPX''), U E AdK{-iHd), (14) 

H 

m = -H - mm (15) 

Let M{t) = P{t)X^(t). The differential equation for M{t) is 

M{t) - [H{t), M(t)], (16) 
such that M{t) E p and the result follows. Q.E.D. 



Remark 7 In the following theorem, we will use the maximum principle, to solve the time optimal 
problem of steering the adjoint control system ( [l0| ) from P(0) = / to the coset KUp, where Up — 
exp{—i9Iizl2zl3z), & [0,47r]. We hasten to add that the proof presented here only establishes that 
the control laws and the corresponding trajectories, given in the following theorem are extremal 
trajectories for the problem of time optimal control. A complete proof of optimality is beyond the 
scope and aim of the present paper and will be presented elsewhere. We first state a lemma which 
will be used in the following theorem. 

Lemma 2 Let A, B, C be as in the notation ||. Then 

exp(27rC) exp(Q;i A + + asC) = I 

for = (2^)'- 
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Proof: First note that exp{tA),exp{tB),exp{tC) are all periodic with period Att and satisfy the 
commutation relation 

[A, B] = C , [B, C] = A, [C, A] = B. 

The mapping A — s- ~ilx, B — s- —ily, C — > —ilz, defines a diffeomorphism between the group 
exp{^, B, C} and SU(2) given by 

exp(aiA + a2B + a^C) exp(-i[Q!i/a; + a2ly + a^Iz]). 

Now using the fact that if X]i=i '^i — (^t'')'^, then eyi]i{—i[ailx +a2ly +ct3lz]) — —I, we obtain that, 
ifEli«? = (27r)2, then 

exp(-j27r/^) exp(-i[ai4 + a2ly + asl^]) = I. 

Therefore 

exp(27rC) exp(aiyl + + a^C) = I. 

Q.E.D. 

Theorem 5 Let Up = expi-iehzhzhz), € [0,47r] and /3 = 27r - 9/2. The control law 

Bt Bt 

H{t) = -i2^J[{hzl2x+hxhz)cOs{^) - ihzl2y+ hyhz) Sm{^)], 

steers the adjoint system ^ from P(0) = / to P{T) e KUp, in 

27rJ 2J 

units of time and is time optimal. 

Proof: Let A,B,C,D be as in the notation ^j. Then using the commutation relations for these 
operators and the BCH, we can rewrite as 

7Y(t) = 27rJexp(-^)^ exp(^). 

The corresponding trajectory P{t), takes the form 

P{t) = exp(-^) exp((^ + 2nJA)t). 

This can be verified by just differentiating the expression for P{t). Next observe that P{T) G KUp- 
To see this note that 

exp(-27rC) exp(27r JTA + /3C) = /, 

where I is the identity matrix. This identity follows directly from the fact (27r JT)^ + = (27r)^ 
and lemma |2[ Therefore 

P(T) = exp(^) = eM-iO{hzhzhz + ^)), 
implying -P(T) G KUp. To see that the control law 7i(i) is extremal, observe for 

M(t) = -7^(t) - ^D, 

the pair {P(t), M{t),'H{t)) satisfies the variation equations |l^, and ^ of theorem |^. To see this, 
recall 

n{t) = 2nJ{Aco,{^)~Bsin{^)), 
12 



therefore the commutation relations 



[A,-D]^B, [B,^D]^-A, 

imply 

[H,M] = ^^[Asin(^) + Scos(|,)]. 

Furthermore 

= ^[Asin(|,) + i3cos(^)]. 

Therefore M{t) satisfies the variational equation M — [7i, M\ and clearly 'H{t) maximizes the func- 
tion tr{HM{t)) for H £ AdKi-i27Tj{IiJ2z + hzhz)) and M{t) = - fz?, Q.E.D. 



Corollary 1 Let [7^? — exp(— i 0/10/2 /3/3-y), G [0,47r] and (a,/?, 7) S {x,y,z). The minimum time 
T, required to steer the adjoint system from P{0) = I to P{T) E KUp, is 



Proof: The proof follows from the observation that hahish-y belongs to the same coset as hzhzhz- 
Therefore the result of theorem ^ apply. 

Proof of Theorem The proof is now a direct consequence of the equivalence theorem ^ and 
theorem ^. 

Geodesic Pulse Sequence: The pulse sequence that produces the propagator 

Up = ejipi-iOhzhzhz), 

in theorem is as follows. 

TT f3 P TT 

Up = cxp{-i-l2y) exp(-i[7r + -]l2x) exp(r(-i27rJ(/i^/2z + hzhz) + ij^hx)) ex:p{i-hy)- 

Where (3 and T are as defined in the above theorem ^. In Fig. |2|(C) a possible implementation of 
this geodesic pulse sequence is schematically shown. Although the simple implementation shown in 
Fig. H(C) is constrained in terms of bandwidth, it forms the basis of more broad band sequence 
which will be presented in a future experimental paper. 



6 Indirect Swap Gates and Coherence Transfer in 3- Spin 

Networks 



In this section, we will consider the problem of transfer of in-phase coherence/^ to /g , for the 
heteronuclear three spin network described by the equation (|^. 

Lemma 3 The unitary propagator 

Vf = exp(-i27r(/i^/2z/3z + hyhzhy + hxhzhx)), 
completely transfers the coherence to . 
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Proof: First observe that Iizhzhz, Iiyl2zhy, and hxhzlsx commute, therefore 

Vf = eyip{-i2iTlizhzhz) exp(-i27r/iy/2z/3y) e^Y>{-i'^T^hxhzhx) ■ 
Furthermore, observe that {/ik, '^Iiyhzhzi '^hzhzhz} forms a so(3) Lie algebra. Therefore, 

e-xp{~i-{AIiJ2zhz)) hx exp(i-(4/i^/2^/3^)) = -^hyhzhz- 
Also note that {^Iiyl2zl?,z, ^^Iiyhzhy, hx} forms a so(3) Lie algebra. Therefore 

TT TT 

exp(--i-(4/iy/2z/3y)) -^hyhzhz exp{i-{4:hyhzhy)) ^ hx- 

Combining the above equalities we obtain VpIixVp = hx- Similarly one can ver ify that VphyV^ = 
hy- Hence the proof. Q.E.D. 

Proof of Theorem ^ (Coherence Transfer:) We need to compute the minimum time required 
to produce the propagator 

Vf = ex.p{-i2Trhzhzhz) e^p{-i2Trhyhzhy) eyi]i{-i2nhxhzhx)- 

We have already shown that the minimum time required to produce a propagator of the form 
exp(-i27r/iQ/2/3/37), where (a,/3,7) G {x,y,z) is 



^27r(27r) ~ (Tr)^ _ ^3 
2^ ~ 2J' 

Therefore Vf can be produced in time less than or equal to (see following remark) . Since there 
might be other unitary propogators, that might achieve this coherence transfer and take less time 
to synthesize, we can only claim that the minimum time required to transfer the coherence /j~ to 
is less than or equal to 

Pulse Sequence: The pulse sequence that produces the propagator 

Vf = exp(-i27r(/i^/2z/3z + hyhzhy + hxhzhx)), 

is as follows. Let Ui = ex.p{~i2Tr{hzhzhz)), U2 = exp(-z27r(/iy/2z/3y)) and C/3 = exp(-i27r(/i^/2^/3a; 
Then 

Ui = exp{-i-hy) exp(-z[7r + -]/2a;) exp{T{-i2TrJ{hzhz+hzhz) + i7^hx}) exp{i-hy)- 

TT TT TT TT 

U2 = exp(i-/ij;) e-xpii-hx) Ui eyip{~i-hx) exp{~i-hx)- 

TT TT TT TT 

U3 = exp(~i-/3j,) exp(-z-/ij^) Ui exp{i-hy) exp{i-hy)- 

Finally 

Vf = Ui U2 U3. 

Where /3 = -tt and T = ^3/2 J. 

Remark 8 It can in fact be shown, that the minimum time required to produce the propogator Vp 
in the above theorem is A rigorous proof is beyond the goals of the present paper, however the 
key observation is that, hzhzhz, hyhzhy ^ and hxhzhx commute, therefore the minimum time 
required to produce the propagator 

Vf = exp{-i27Thzhzhz) exp>{-i2TT hyhzhy) cx.p{-i2'n-hxhzhx), 

is the sum of minimum time required to produce the individual propagators exp(— i2n hzhzhz) j 
eyip{-i2TThyhzhy) and exp(-i27r/i3;/22-^32;). 
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Proof of Theorem |^: (Indirect Swap Gates) The indirect swap gate Usw{^, 3) is given by 

TT 

Usw{l, 3) = exp{-i2n{Iizl2zl3z + hyhzhy + hxhzhx)) e-xpii-hz)- 

The propagator exp(i^/2z) can be produced in arbitrarily small time by selective hard pulses. 
Therefore the minimum time required to produce the swap gate is the same as the minimum time 
required for creating exp(— i27r(/iz/2z/3z + Iiyhzhy + Iixhzhx)), which is Hence the proof 

Q.E.D. 



Remark 9 Synthesis of A2{U) gates: Pulse sequences for produce A2 gates, in the context of 
NMR quantum computing need to synthesize effective Hamiltonians of the form /iq/2/3^37- To see 
this, observe that 

" 1 
1 



A2(4 







1 

1 



1 

1 









1 












-1 



This can be rewritten as 



A2(4) = exp(-i7r 



Thus the effective Hamiltonian takes the form 



" 


■ 




■ 


■ 




■ 


" 





1 







1 







1 



= exp[-i7r(y - Ji^) «> (y - hz) ® (y ^ -^32) 



h\ - hz) ® ( Y - hz) 



7r(. 



hz + 1: 



3z) 



(y-/3.) 

{hzhz + hzhz 



Iizhz) 
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hzhzhz ) 



Since the term IizhzHz commutes with other terms in the effective Hamiltonian, it needs to be 
produced besides the other terms in the -ffe// to synthesize the A2 (Iz) gate. We have already com- 
puted the time optimal pulse sequences for the optimal implementation of an effective Hamiltonian 
of the form Iizhzlsz- Therefore to derive optimal implementations of A2(/z) gates, further work 
is required to compute is shortest pulse sequences for synthesizing an effective Hamiltonian of the 
form /iz/sz- 



7 Conclusion 



In this paper we have demonstrated substantial improvement in the time that is required to syn- 
thesize an important class of unitary transformations in spin systems consisting of three spins ^. It 
was shown that computing the time-optimal way to transfer coherence in a coupled spin network 
can be reduced to problems of computing sub-Riemannian geodesies [^4j . These problems were then 
explicitly solved for a linear three spin chain. These ideas are not just restricted to the 3-spin case 
considered in this paper but can be extended to find time optimal pulse sequences in a general 
quantum network [ p^ . 
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Table 



Comparison 



of 



Pulse 



Sequence 



Durations 



Unitary Transformation T(State of the art sequences) r* (Geodesic sequences) 



Up = exp{-i27rfilial2pl3'r) 
Up = exp(-i27r/iQ/2/3/37) 
Swap(l,3) 



2 + K 

2J 
_3_ 
2J 
_9_ 
2J 

3. 

■ I 



Vl>^l(4-|>c|) 
2\J\ 
V3 

2171 

3\/3 
2J 
3%/3 
21 



\/I«^I(4~|k|) 

2+fc 

73 = 5^-^^° 
^ = 86.6% 



A Appendix: Spin-decoupling 



Given the evolution of the unitary propogator 

m 

IJ = -i{Hd + J2 ^iHj) U, U{0) = I, 
i=i 

let have a decomposition — Hf + such that [Hf^H^] — 0. The control Hamiltonians 
{Hj}, generate a subgroup K, given by 

K = exp({i/-,}LA), 

where {Hj}LA is the Lie algebra generated by {—iHi, — ii?2, ■ • • , —iH„i}. Let fc G A' be such that 

fc-i - ^(i/,^ + i7|)fc = -.(i/,^ - i/f ). (17) 

It is assumed that the strength of the control Hamiltonians can be made arbitrary large. Under this 
assumption the propogator k can be produced in arbitrarily small time, such that the evolution due 
to the drift Hd during this time can be neglected. Now consider the evolution 

U{t) = exp(-ii/d^) fc"^ exp(-ii/<j^) k. 

From equation |l^, we obtain 

U{t) = cxp{-z[H^ + Hi]i) cxp(-z[i7f - Hi]i) = cxp(-zi7f i)- 

Therefore the net evolution is as if the system evolved under the drift term for time t. We 
will say that the part of the drift has been decoupled. In a network of coupled spins, 
represents the coupling of a specified spin to the rest of the network and decoupling corresponds 
to decoupling the spin from the network. 
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Figure 2: The panel shows three pulse sequences for synthesizing the trilinear propagator Up = 
exp{—i9Iizl2zl3z) with 9 = 2'kk. The conventional pulse sequence A uses decoupling and takes time 
t = (2 + k)/2J. The second pulse sequence B improves the first sequence by avoiding decoupling 
and has a duration <' = (! + K)/2J.The final pulse sequence C is time optimal and has a duration 
t* = \/k(4 — k)/2J. The radio-frequency amplitude Vrf of the hatched pulse is (2 — /t) J/-\/k(4 — k). 
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